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ABSTRACT 

We report new mid-eclipse times of the short-period sdB/dM binary HWVirginis, which differ substantially from the times predicted 
by a previous model. The proposed orbits of the two planets in that model are found to be unstable. We present a new secularly stable 
solution, which involves two companions orbiting HW Vir with periods of 12.7 yr and 55 ± 15 yr. For orbits coplanar with the binary, 
the inner companion is a giant planet with mass Ms sin (3 =i 14 Mj up and the outer one a brown dwarf or low-mass star with a mass of 
M4 sin (4 = 30 - 120 Mj up . Using the mercury6 code, we find that such a system would be stable over more than 10 7 yr, in spite of 
the sizeable interaction. Our model fits the observed eclipse-time variations by the light-travel time effect alone, without invoking any 
additional process, and provides support for the planetary hypothesis of the eclipse-time variations in close binaries. The signature of 
non-Keplerian orbits may be visible in the data. 
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1. Introduction 

Periodic or quasiperiodic variations of the mid-eclipse times in 
close binaries have been observed for decades and variously as- 
cribed to activity cycles of the secondary star (Applegate 1 992), 
apsidal motion (Todoran 1972), or the response to a third body 
orbiting the binary (Nather & Robinson 1974). More recently, 
Lee et al. (2009) assigned the complex O-C (observed minus 
calculated) eclipse-time variations in the detached sdB/dM bi- 
nary HW Vir with a 2.8-h orbital period to the light-travel time 
(LTT) effect caused by the orbital motion of two giant plan- 
ets. Thereafter, an increasing number of eclipsing post-common 
envelope binaries (PCEB), both with sdB and white-dwarf pri- 
maries, were proposed as having eclipse-time variations pos- 
sibly due to planetary or brown-dwarf companions. Systems 
with one proposed companion include HS0705+67, DP Leo, 
HS2231+24, and NSVS 14256825 (Qian et al. 2009, 2010a,b; 
Beuermann et al. 2011, 2012). Further PCEB that may harbor 
more than one companion are NN Ser (Beuermann et al. 2010), 
UZ For (Potter et al. 2011), RRCae and NY Vir (Qian et al. 
2012a,b), QS Vir (Parsons et al. 2010), and HU Aqr (Qian et al. 
2011; Hinseetal. 2012). 

The question of secular stability does not arise for binaries 
with a single planet, but must be considered for the proposed sys- 
tems of companions. Such analyses were not included in most 
original publications, with the exception of the two planets that 
orbit NN Ser, which represent a stable resonant pair (Beuermann 
et al. 2010). The system suggested by Qian et al. (2011) for 
HU Aqr, on the other hand, was shown to be secularly unstable 
(Horner et al. 2011; Wittenmyer et al. 2011), although a differ- 
ent stable solution may exist (Hinse et al. 2012). To be sure, the 
planetary hypothesis of the observed eclipse time variations has 
met with some scepticism (e.g. Wittenmyer et al. 2011), but a 
sufficiently well-defined alternative is not in sight. 

A more general critique of the models of Lee et al. (2009) 
and Qian et al. (2011) relates to their combining the LTT ef- 



fect with a long-term period decrease on a time scale of t = 
P/P o£ 10 7 yr, supposedly produced by some other mechanism. 
Gravitational radiation leads to a period decrease, but is entirely 
negligible with r 3 x 10 9 yr. Magnetic braking is a more ade- 
quate contender for a long-term period decrease, but cannot ex- 
plain the recent period increase in HW Vir. Applegate's (1992) 
mechanism allows variations of both signs, but most authors 
agree that it is too feeble to produce the large observed eclipse- 
time variations (e.g. Brinkworth et al. 2006; Watson & Marsh 
2010). If, however, a so far unknown super-Applegate effect 
were operative in close binaries, it could conceivably account 
for the entire observed effect and render the arbitrary division 
into variations produced by different mechanisms meaningless. 
Indeed, it is the spectre of a super-Applegate mechanism that 
has historically prevented us from considering in close binary 
systems what, for unevolved single stars, has become the default 
situation: complex planetary systems (Lovis et al. 2011). The 
discovery by KEPLER of several circumbinary planetary sys- 
tems in non-evolved eclipsing binaries (Doyle et al. 201 1 ; Welsh 
et al. 2012) shows us that it is not unreasonable to consider the 
same for evolved binaries, although their planetary systems pos- 
sibly do not not survive the common-envelope phase unscathed. 
In this light, the proposed circumbinary planetary systems de- 
tected by the LTT method must be carefully re-scrutinized. 

Here, we present new mid-eclipse times for HW Vir that de- 
viate significantly from the Lee et al. (2009) prediction. We show 
that the planetary system proposed by Lee et al. is secularly un- 
stable and, therefore, untenable. We find that a secularly stable 
two-companion model can be devised, in which the observed 
eclipse-time variations are due to the LTT effect alone. 



2. The data base 

Starting in December 2010, we monitored the V=10.9 mag bi- 
nary HW Vir with the MONET/North telescope at the University 
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Table 1. New primary mid-eclipse times T xl of HWVir and 0-Cu n> i 
residuals for the MONET/N data relative to the ephemeris of Eq. 1 . 



Cycle 


BJD(TT) 


Error 


O— C]j n> i 


Origin 




2400000+ 


(days) 


(days) 




75123 


54498.877674 


0.000014 


N.A. 


AAVSO 


75890 


54588.401364 


0.000040 


N.A. 


AAVSO 


76004 


54601.707367 


0.000016 


N.A. 


AAVSO 


76063 


54608.593786 


0.000019 


N.A. 


AAVSO 


76089 


54611.628553 


0.000012 


N.A. 


AAVSO 


78062 


54841.916149 


0.000010 


N.A. 


AAVSO 


84077 


55543.984048 


0.000008 


-0.000007 


MONET/N 


84120 


55549.003005 


0.000008 


0.000009 


MONET/N 


84180 


55556.006176 


0.000009 


0.000006 


MONET/N 


84411 


55582.968393 


0.000009 


0.000006 


MONET/N 


84428 


55584.952622 


0.000010 


0.000003 


MONET/N 


84488 


55591.955807 


0.000010 


0.000014 


MONET/N 


84497 


55593.006274 


0.000007 


0.000005 


MONET/N 


84600 


55605.028372 


0.000007 


-0.000010 


MONET/N 


84608 


55605.962117 


0.000015 


-0.000022 


MONET/N 


84693 


55615.883298 


0.000007 


-0.000004 


MONET/N 


84863 


55635.725619 


0.000006 


-0.000007 


MONET/N 


Q A OA"7 


jju4/.ou44oU 


U.UUUUU/ 


U.UUUUU 1 


JYLUJNrll/JN 


84976 


55648.914932 


0.000007 


-0.000004 


MONET/N 


85026 


55654.750921 


0.000006 


0.000007 


MONET/N 


85249 


55680.779371 


0.000007 


-0.000003 


MONET/N 


85266 


55682.763597 


0.000015 


-0.000009 


MONET/N 


87093 


55896.010239 


0.000007 


0.000005 


MONET/N 


87589 


55953.903110 


0.000018 


-0.000023 


MONET/N 


87624 


55957.988315 


0.000007 


-0.000003 


MONET/N 


87787 


55977.013609 


0.000007 


0.000004 


MONET/N 



of Texas' McDonald Observatory via the MONET browser- 
based remote-observing interface. The photometric data were 
taken with an Apogee ALTA E47+ lkxlk CCD camera mostly 
in the / c -band with exposure times of 10 or 20 s. Since there is 
no suitable comparison star in our 5' x 5' field, we included only 
observations obtained under photometric conditions. The light 
curves were analyzed using the heuristic mathematical model 
described in paper II of this series (Beuermann et al. 2012). 
Table 1 lists the 20 new primary mid-eclipse times along with 
their formal 1 -<x errors, which range from 0.5 to 1 .6 s. Fitting the 
MONET/N mid-eclipse times alone yields the linear ephemeris 
valid in 2010/2012 

r ed = BJD(TT) 2455543.984055(2) + 0.1 16719555(2) E. (1) 

The residuals from Eq. 1 are listed in Table 1 as 0—C\i n ^ . Their 
rms value of 0.8 s is consistent with the errors obtained from the 
formal fits to the light curves. Because of the effects of the addi- 
tional bodies in the system, the period of Eq. 1 is not necessarily 
identical to the binary period. 

A large body of primary and secondary mid-eclipse times 
of HW Vir is available in the literature. The S AAO group mon- 
itored HWVir between 1984 and 2002 (Marang & Kilkenny 
1989) and Kilkenny et al. (1991, 1994, 2000, 2003), reporting 
a total of 1 1 1 primary mid-eclipse times with errors mostly as 
small as 2 s. The measurements of Lee et al. (2009) overlap with 
the SAAO data and extend the coverage to 2009, with a gap in 
2004. Additional mid-eclipse times were published by Wood et 
al. (1993), Giirol & Selan (1994), Wood & Saffer (1999), Cakirli 
& Devlen (1999), Kiss et al. (2000), Agerer & Hiibscher (2000, 
2002, 2003), ibanoglu et al. (2004), and Brat et al. (2008, 2009, 
2011). Further timings were drawn from the Japanese VSNET 
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Fig. l.O — Cii n .2 residuals of the mid-eclipse times from the linear 
ephemeris used by Lee et al. (2009) along with their model curves 
for the two-companion model (solid) and the underlying quadratic 
ephemeris (dashed). The data are from SAAO (green), Wood et al (cyan 
blue), Lee et al. (yellow), BAV and VSNET (magenta), AAVSO (blue), 
BRNO (red), and MONET/North (green). 



archive (Kato et al. 2004) and unanalyzed light curves were 
obtained from the AAVSO archive (Henden 20 10) 2 . We deter- 
mined mid-eclipse times for the AAVSO data, using the same 
fitting method as for the MONET data. The resulting new pri- 
mary mid-eclipse times are included in Table 1 . 

The general picture that emerges from the entire body 
of eclipse times collected between 1984 and 2012 is that of 
a smooth long-term O - C variation of complex shape (e.g. 
Kilkenny et al. 2003; ibanoglu et al. 2004; Lee et al. 2009). Our 
own observations demonstrate the absence of short-term O — C 
variations with periods less than 1 yr exceeding a couple of sec- 
onds. We have scrutinizingly surveyed the available data and 
found that a small number of published mid-eclipse times de- 
viate significantly from the mean O—C variation, suggesting that 
the errors were underestimated. Rather than including all data 
indiscriminately, we excluded these outliers from our analysis 
on the assumption that the absence of short-term O—C variations 
as documented by the SAAO timings and our own data holds 
at all times. In general, we accepted mid-eclipse times if they 
have a quoted error not exceeding 0.0001 day. An adopted sys- 
tematic error of 3.0 s was quadratically added to all eclipse times 
included in our data base. We restricted our analysis to primary 
mid-eclipse times, because the times for the secondary eclipses 
are significantly less well-determined and add little to the defi- 
nition of the long-term O-C variation. The mean phase of the 
secondary eclipses agrees with <p - 0.500 within a few seconds, 
similar to our finding for the sdB systems NSVS 14256825 and 
HS 0705+67 (paper II). 

Figure 1 shows the 0-C\ m ,7. residuals of the data set adopted 
by us relative to the linear ephemeris used by Lee et al. (2009) in 
their Fig. 5 (top panel). We have included all 111 primary mid- 
eclipse times reported by the SAAO group, 20 from Lee et al. 
(2009), four from Wood et al. (1993) and Wood & Saffer (1999), 
three from the OEJV (Brat et al. 2011), eight VSNET times, 
three BAV times, and the new mid-eclipse times of Table 1. 
Our entire data base contains 176 primary mid-eclipse times. 
We have corrected all times to the Solar-system barycenter in 
the terrestrial system quoting them as BJD(TT) 3 . 



1 http://www.kusastro.kyoto-u.ac.jp/vsnet/ 

2 http://www.aavso.org/ql 

3 http://astroutils.astronomy.ohio-state.edu/time/ 
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Table 2. Parameters of the eclipsing binary HW Vir. 



Parameter 



Value ± Error 



Reference 



Orbital period _P bin (s) 
Primary mass M x (M Q ) 
Secondary mass M 2 (M Q ) 
Secondary radius R 2 (R B ) 
Separation a bin (Re) 
Inclination ; bin (°) 
Eccentricity ebin 
Distance d (pc) 
Visual magnitude V (mag) 



10084.58 +0.01 
0.485 ±0.013 
0.142 ±0.004 
0.175 ±0.026 
0.860 ±0.010 
80.9 ±0.1 
< 0.0003 
181 ±20 
10.9 



This work 
Lee et al. (2009) 
Lee et al. (2009) 
Lee et al. (2009) 
Lee et al. (2009) 
Lee et al. (2009) 
This work 
Lee et al. (2009) 




35 

period (yr) 

Fig. 2. Contour plot of the reduced x 1 normalized to unity for the best 
fit using a generalized Lomb-Scargle periodogram for a two-companion 
fit to the eclipse-time variations of HW Vir. The solid lines indicate 2:1, 
3: 1 and 4: 1 resonant orbits, the cross the best fit. Solutions with an inner 
planet of about 12.7 year orbital period provide good fits for a wide 
range of orbital periods of the outer companion. 



In Table 2, we summarize the parameters of the binary 
HW Vir relevant to the present study. The masses and the dis- 
tance are taken from Lee et al. (2009). The limit on the ec- 
centricity was obtained from our new mid-eclipse times r ec i in 
Table 1, which limit the amplitude caused by apsidal motion 
to Ar ec i =s P hin e/n < 1.0 s for a period of the apsidal rotation 
U * P bi n(M 2 /M l )(a bin /R 2 ) 5 /(l5k2) * 43d, with k 2 * 0.15 the 
structure constant of the nearly fully convective secondary star 
(Feidenetal. 2011). 

3. The Lee et al. model 

In an influential paper, Lee et al. (2009) interpreted the mid- 
eclipse times of HWVir available until 2008 (cycle number 
E = 76050) by the LTT effect of two planets superposed on a 
quadratic variation of unspecified origin. Figure 1 shows their 
model curves for the quadratic variation (dashed curve) and their 
final model (solid curve). Their fit is adequate until 2008, but 
fails completely to reproduce our new data. The discrepancy has 
reached 250 s in 2012 and is largely due to the rapid fall-off of 
the quadratic term. The continuous period decrease as defined 
by the quadratic term does not exist. 

The Lee et al. (2009) model faces the additional problem that 
the proposed 2-planet system is secularly unstable. This result is 
evident from the planetary parameters quoted by them, which 
imply crossing orbits with an apoapsis of the inner planet of 
4.7 AU and a periapsis of the outer planet of of 2.9 AU. We have 
numerically integrated the orbits with mercury6 (Chambers 
1999, see below for more details) and find that a near encounter 
or a collision occurs within 2000 yr. Hence, the model of Lee et 
al. (2009) is untenable in the present form. 

In the remainder of the paper, we show that all observations 
can be explained by the LTT effect, without taking recourse to 
an additional mechanism. A minor contribution by gravitational 
radiation is not excluded, but is below our detection limit. 

4. Search method 

We first scanned the two-companion parameter space for orbital 
periods below 60 years, using a generalized Lomb-Scargle peri- 
odogram supplemented by refined local searches. In a second 
step, we searched for improved solutions using a Levenberg- 
Marquardt (LM) minimization algorithm allowing us to de- 
tect the true local minima near the start parameters. Finally, we 
tested the secular stability of all possible solutions using the hy- 
brid symplectic algorithm implemented in the mercury6 pack- 
age (Chambers 1999). 



Fitting Keplerian orbits to the set of mid-eclipse times in- 
volves a simultaneous fit for the ephemeris of the binary star 
(two parameters) and the orbits of the companions (five pa- 
rameters each). The orbital inclination i of a companion re- 
mains undetermined. For two or more companions, a direct 
search of the parameter space is computationally elaborate and 
time-consuming because of the large number of parameters in- 
volved. We adopted, therefore, the approach of Zechmeister & 
Kiirster (2009) and developed a variant of the generalized Lomb- 
Scargle (GLS) periodogram (Lomb 1976; Scargle 1982; Home 
& Baliunas 1986). The mid-eclipse time of eclipse number n in- 
cluding the light-travel time of planets k = 1 . . . N p is expressed 
as 

l-e 2 

T(n) = T + n P bin + V K k k —— sin (v k (n) + oj k ) , (2) 

f^f l+e k cosv k (n) 

where e k is the eccentricity of planet k, K k = a b i ak sin i k /c the 
amplitude of the LTT effect, with a b \^ k the semi-major axis of 
the orbit of the center of gravity of the binary around the com- 
mon center of mass and i k the inclination, oj k is the argument of 
periastron measured from the ascending node in the plane of the 
sky, and v k ( n ) the true anomaly at time t — To + n Pbin (Kopal 
1959). Eq. 2 can be transformed into a linear equation for the 
coefficients A k , B k , and C, using the eccentric anomaly E k (n) of 
planet k at the time of eclipse n, 

T(n) = n Pbi„ + C + J^(A k cos E k (n) + B k sin E k (nj) (3) 

k=\ 

with 



A k = -K k sinaj k , B k = K k ^l-e z k costi> k , C = r +2 J K k e k sinw^. 

k=\ 

The best-fit values of Tq, P b { n , K k , and u k were derived by linear 
regression for given values of e k , orbital period P k and time t k of 
periastron passage. Calculating a grid in these quantities reduces 
the number of free parameters for each grid point from 2 +5N V 
to 3 N p . 

In a second step, we searched for improved solutions near 
the best solution from GLS scan. We employed the Levenberg- 
Marquardt (LM) minimization algorithm (Markwardt 2009), a 
non-linear least-squares fitting routine implemented in mpfit 
of IDL. We used the GLS-parameters as start values and gener- 
ated 500 variations of them, using the diagonal elements of the 
mpf it-derived covariance matrix as uncertainties cr. The start 
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Cycles 

Fig. 3. Fit of two Keplerian orbits to the eclipse-time variations of 
HW Vir. Top: Data of Fig. 1 relative to the linear ephemeris of Eq. 4. 
The curves denote the model LTT effect (solid) and the contributions by 
the outer companion (long dashes) and the inner planet (short dashes). 
Center: Data with the contribution by the outer companion subtracted 
and model for the inner planet (solid curve). Bottom: Residuals after the 
subtraction of the contributions by both companions. 

parameters were varied randomly within +10 cr around the orig- 
inal start values, allowing us to detect the true local minimum 
near the start parameters. 

Finally, we tested the secular stability of the solutions, using 
the hybrid symplectic integrator in mercury6 (Chambers 1999), 
which allows one to evolve planetary systems with high pre- 
cision over long times very efficiently. We used constant time 
steps of 35 d. Test runs demonstrated that this choice is adequate 
and smaller steps do not change the results. This provision is 
not adequate for the treatment of a close encounter, but such 
incidence should not occur in the successful models and if it 
does, the calculation is stopped and the model termed 'unsta- 
ble'. The accuracy parameter in the code was set to 10~ 16 , lead- 
ing to a fractional change in the energy and momentum of the 
triple (binary and two companions) of typically 10~ 5 - 10~ 4 and 
10~ 9 - 10~ 8 , respectively. The quoted changes are valid for an 
integration time of 10 7 yr and the larger values obtain for longer 
periods P4 and, hence, larger masses M4. For simplicity, the cen- 
tral binary was treated as a single object with a mass equal to the 
sum of the component masses. This simplification is justified 
given the short binary orbital period of 2.8 h. The gravitational 
field at the position of a distant companion can be represented as 
the sum of the constant field created by the combined mass of the 
binary components and a gravitational wave, emanating from the 
revolving binary with periods of 2.8 h for the fundamental and 
1.4h for the first harmonic. The relative strength of the wave 
field is 4 X 10~ 8 and closely averages to zero over the 300 or 600 
periods that occur in a time step of 35 d. The retroaction of the 
companion tends to excite an eccentricity in the binary, but with 



a relative strength of 10 this effect is also entirely negligible. 
We used the masses of the binary components as given by Lee 
et al. (2009) and quoted in Table 2. The orbital evolution of the 
companions was followed until instability occurred, or at least 
for 10 7 yr and up to 10 8 yr for some models. 

5. A stable two-companion model for HW Vir 

The data shown in Fig. 1 suggest the presence of a period near 
40000 cycles or 13yr superposed on a variation with a longer 
period. Using the GLS periodogram, we scanned the parameter 
space for two companions with orbital periods of P3 = 10- 15 yr 
and P4/P3 > 1.85. Figure 2 shows a contour plot of the reduced 
X 2 normalized to unity for the best fit. For a better presentation, 
only the range P4 <45 yr is displayed; no new features appear for 
larger P4. A ridge of low x 2 confirms P3 ^ 12.5 - 12.9 yr for all 
P4 ^ 28 yr, including the possible mean-motion resonances with 
P4 : P 3 = 5 : 2, 3 : 1, 7 : 2, 4 : 1, and 5 : 1. None of the resonant 
solutions, however, is preferred over non-resonant ones. For the 
best fit near P4 = 30 yr , the inner object has a mass close to the 
boundary between planets and brown dwarfs, the outer one is a 
brown dwarf, with a mass increasing with orbital period up to 
the stellar mass limit for P4 70 yr. 

With the LM optimization algorithm, we found improved 
two-companion solutions along the ridge of low x 2 m Fig- 2. 
Formally, the best fit is now obtained near a period ratio of 4 : 1 
with periods P3 - 12.7 yr and P4 around 50 yr. The fit deterio- 
rates for P4 > 70 yr. The eccentricity of the inner companion is 
rather well defined with =0.4+0.1, that of the outer one is not, 
although small finite values fit better than e4 = 0. Since the data 
cover only part of the orbital period of the outer companion, the 
deduced amplitude K4 of the LTT effect depends on the adopted 
values of P4 and e/j. Aided by the stability calculations described 
below, we chose of P4 = 55 yr and e4 = 0.05. Figure 3 (top panel) 
shows the best fit and Table 3 summarizes the fit results with our 
best estimates of the errors. The center panel repeats the data 
with the dominant contribution from the outer companion sub- 
tracted (<9-C e ii,i) and the bottom panel shows the residuals after 
subtraction of both contributions (O - C e ii,2)- Small systematic 
residuals exist, notably around cycle number 70000, where the 
data quality is not the best, however. The formal best fit leads 
to a symmetric configuration with Ao> = 0)4-0)3 near zero, but 
with a large uncertainty in 0)4, because of the limited data cov- 
erage. The general decrease of 0-Cn n over the 1984-2012 time 
interval reflects the fact that the observations cover only that part 



Table 3. Parameters of the two-companion model of Fig. 3 for HW Vir. 
A colon indicates uncertain values. 



Parameter 


Inner companion 


Outer companion 




#3 




#4 


Orbital period P (yr) 


12.7 ± 


0.2 


55 (fixed) 


Eccentricity e 


0.40 ± 


0.10 


0.05 : 


Semi-major axis a (AU) 


4.69 ± 


0.06 


12.8 + 0.2 


Amplitude K (s) 


49 ± 


3 


563 ± 200 


Mass Msin; (M Jllp ) 


14.3 ± 


1.0 


65 ± 15 


Argument of periastron a> ( c 


) -18 + 


10 


0: 


Periastron passage (JD) 


2452401 


2461 677 : 


Periastron passage (Cycle) 


57 150 


136619 : 



4 



20 




° 10 

G 

E 5 

CD 

n 

0.5 

* 0.4 

CO 
CD 

£ 0.3 



0.2 




0.0 
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Fig. 4. Temporal variation of the semi-major axes and the eccentricities of the osculating orbits of the companions to HW Vir as calculated with 
the Mercury6 code (see text). The system is stable for more than 10 8 yr (see text). Color coding is the same for all panels. 



of the orbit, when the outer companion is receding from the ob- 
server and the binary is approaching. The corresponding LTT 
amplitude suggests that the outer companion is probably a brown 
dwarf with a mass of roughly 65 M Jup or a star of very low mass. 
For periods between 40 and 70 yr, its mass would be in the range 
of 30 to 120Mj up . The fit yields an underlying linear ephemeris 
of the binary 

r ed = BJD(TT) 2445730.5497(72) + 0.1 1671969(15)£, (4) 

where the quoted correlated uncertainties of the epoch and the 
binary period Pbin refer to the quoted range of P4. 

Using the mercury6 code, we followed the orbital evolution 
of solutions along the low-^ 2 ridge in Fig. 2 numerically, us- 
ing combinations of periods P3 and P4 and masses M3 and M4 
that match the observations. We note that the actual positions of 
the two companions in space, their orbital velocities, and true 
anomalies at a given time can not be obtained with sufficient 
accuracy from the observations. We employed, therefore, stan- 
dardized start parameters and searched for stable models that 
reproduce the observations for a substantial percentage of the 
time. The calculations were started in either the symmetric or 
the antisymmetric configuration with both companions at peri- 
astron, having Keplerian velocities. The orbits were taken to be 
coplanar with the binary orbit (2 = 80.9°, Table 2), with the start 
eccentricities as free parameters. Orbital evolution changes the 
elements quickly and periodically varying eccentricities develop 
also for models that start with circular orbits. 

We dubbed a model 'stable' if it persisted without close en- 
counter for more that 10 7 yr. Such models are found for periods 
P4 as short as 35 yr, but for most start configurations at this pe- 
riod instability incurred, sometimes only after 10 6 yr. The frac- 
tion of stable models that mimic the observations becomes sub- 
stantial only near 45 yr and at 55 yr or longer. A low fraction 
of stable models near 37 yr and, less so, near 50 yr may imply 
a higher probability of instability near the 3:1 and 4:1 mean- 
motion resonances. The general preference for longer periods 



P4 is expected, because the observations fix Pj between 12.5 and 
12.9 yr and a more distant outer companion increases the phase 
space for solutions with larger eccentricities. At P4 > 55 yr, sta- 
ble models were found for eccentricities <?3 as large as 0.45. The 
properties of stable models between P4 = 45 yr and 55 yr differ 
in that the former prefer anti-symmetric configurations, whereas 
the latter all spend a substantial fraction of the time near a sym- 
metric configuration. Taken at face value, our observational best 
fit of Alo near zero favors a model with P4 >50 yr. Our preferred 
model has P4 = 55 yr, values of P3, e^, M3, and M4 as obtained 
from the fit to the data (Fig. 3, Table 3), and orbits coplanar with 
the binary. 

Figure 4 shows the evolution of the semi-major axes and 
eccentricities for the preferred model. The left panels illus- 
trate the behavior for the first 165 yr (three orbital periods of 
the outer companion), the right panels show the evolution over 
lOOOOyr. The system is stable for 10 s yr. The elements of the 
osculating orbits vary periodically on time scales up to 4000 yr 
(right panels). The eccentricities and 64 range from 0.32 to 
0.43 and from zero to 0.18, respectively. Over short time inter- 
vals, the modulation at the synodic period of the inner planet 
P syn = (Pj 1 - P 4 l y l - 16 yr is dominant (left panels). The sys- 
tem spends a substantial fraction of the time near the observed 
eccentricities, ej, 0.4 and e4 0.05, and assumes a symmetric 
configuration most frequently near maximum in the 4000 yr 
cycle (lower right panel), suggesting that the observational fit of 
Fig. 3 and Table 3 represents a snapshot of the variable consti- 
tution of the HW Vir system. While the agreement seems satis- 
factory, we add the caveat that our simulations have so far only 
provided a first glimpse at the stability landscape, which remains 
to be explored in more detail. 

The perturbation of the orbits by the mutual interaction of the 
two companions leads to non-Keplerian shifts of the observed 
eclipse times, which are not included in the fit shown in Fig. 3. 
Our simulations show that such shifts are not negligible and af- 
fect the derived values of the LTT amplitude, the period, the ec- 
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centricity, and the argument of periastron of both companions. In 
particular, the LTT amplitude of the inner companion can vary 
by more than 10 s between subsequent orbits. Such orbit-to-orbit 
variations are not correctly interpreted by our present fitting rou- 
tines. The observed residuals 0-C e ii,2 after subtraction of the 
contributions from the two Keplerian contributions may show 
a non-random structure, possibly representing the signature of 
the perturbed orbits (Fig. 3, bottom panel), but the quality of the 
data is lowest between cycle numbers 60000 and 80000, where 
the excursion appears largest. Clearly, a unique identification of 
non-Keplerian eclipse-time variations in PCEB would provide a 
strong support for the planetary model. 

We have alternatively considered that the residuals in Fig. 3 
might indicate the presence of a third inner planet. While the 
added free parameters improve the fit, this possibility can be ex- 
cluded, because such object is quickly expelled. 

6. Conclusion 

We have investigated the plausibility of the planetary model for 
HW Vir, using additional data and an exhaustive search for and 
dynamical tests of two-companion solutions. The original model 
by Lee et al. (2009) is clearly ruled out, both from the current 
direction of the observed O - C variations as well as the fact 
that their solution is dynamically unstable. We have presented 
a qualitatively different solution that involves two companions 
in secularly stable orbits about the binary HWVir. The inner 
companion has an orbital period P3 12.7 + 0.2 yr and a mass 
M3 14.3 + 1.0Mj U p for an inclination z'3 identical to that of the 
binary (z = 80.9°, see Table 2). For other inclinations, the mass 
varies as 1/sin z'3. The outer more massive companion is a brown 
dwarf or a low-mass star with a mass around M4 sin z'4 ^ 30 to 
120Mj up , where z'4 is the unknown inclination of its orbit. Its pa- 
rameters are still uncertain, because the data cover only about 
half an orbital period. While the suggested pair of companions 
in coplanar orbits is secularly stable, we caution that the stabil- 
ity landscape has not been thoroughly explored and its system- 
atic structure is still eluding us. The system is strongly interact- 
ing with evidence for non-Keplerian orbits, whose signature may 
even be visible in the data. 

Our model of HW Vir provides an attractive possibility of ex- 
plaining the entire observed orbital-period variations of HW Vir 
by the LTT effect, without taking recourse to any additional un- 
explained process. It provides positive evidence of the planetary 
model for HW Vir and strengthens the case for close binaries in 
general. It also calls for a re-examination of other systems in 
the light of our result. However, the case of HW Vir also shows 
how difficult it may be to find such a solution, given both the 
very long time-series required to detect clearly periodic plane- 
tary signals, the requirement of fitting all residuals without in- 
voking other arbitrary mechanisms, and the limits of fitting in- 
dividual static Keplerian orbits to what must obviously be very 
dynamic systems. 
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